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1.  INTRODUCTION 

This  report  presents  the  results  of  an  investigation  into  the  feasibility  of  using  Disturbance 
Accommodating  Controller  (DAC)  design  techniques,  as  developed  by  Dr.  C.  D.  Johnson  of 
the  University  of  Alabama  in  Huntsville,  to  cancel  out  disturbance  inputs  to  a  missile 
autopilot  channel. 

The  DAC  method  of  design  uses  a  combination  of  waveform-mode  disturbance  modeling 
and  state-variable  control  techniques.  As  a  tool  for  controller  design,  the  DAC  approach 
permits  three  primary  modes  of  disturbance  accommodation:  (I)  cancellation  (absorption)  of 
disturbance  effects,  (2)  minimization  of  disturbance  effects,  or  (3)  constructive  utilization  of 
the  disturbances  as  an  aid  in  accomplishing  the  primary  control  task. 

The  purpose  of  this  report  is  to  determine  if  these  techniques,  specifically  the  cancellation 
and  minimization  modes,  can  be  successfully  applied  to  a  missile  system  in  such  a  manner  as  to 
cancel  out  the  effects  of  disturbance  inputs  which  would  otherwise  degrade  system  accuracy. 

It  is  not  the  intent  of  this  report  to  thoroughly  cover  all  the  background  theory  involved  in 
the  development  of  DAC  design  procedures.  This  theory  can  best  be  obtained  by  reading  the 
original  papers;  see,  for  instance,  References  1-4.  For  applications  of  DAC  to  several  simple 
systems  see  Reference  5. 

2.  SOME  BACKGROUND 

The  plant  considered  in  this  report  is  one  which  can  be  described  by  state  equations  of  the 
form 


x  «  Ax  +  Bu  +  PW 

y  -  Cx  +  Eu  +  Ow_  (1) 


where 

i  is  the  plant  state  vector, 

jl  is  the  plant  control  input  vector, 

w  is  the  vector  of  external  disturbance  acting  on  the  plant, 

y  is  the  system  output  vector,  and 

A.  &  F,  C,  E,  £  are  appropriate  size,  known  matrices  which  are  not  necessarily  constant. 
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Now,  the  external  disturbances,  w(t),  for  which  DA C  theory  is  intended  are  characterized 
by  the  presence  of  “waveform  structure.”  i.e.,  the  functions  w(t)  can  be  described  by  known 
differential  equations  which  the  w(t)  satisfy  “almost  everywhere.”  For  the  cases  considered  in 
this  report,  the  disturbances  will  be  assumed  to  be  described  by  the  following  general  set  of 
linear  disturbance  state  equations: 

H  -  HS.  +XJL 

z  +  De  +  Mx  +  a  (2) 


where 

z_  is  the  disturbance  “state”  vector, 

a  is  a  sequence  of  randomly  arriving  vector  impulses,  and 
J),  H,  L,  M  are  known,  time-invariant  matrices. 


In  most  practical  applications,  neither  the  complete  set  of  plant  state  variables  nor  the 
various  components  wj|(t)  of  the  disturbance  are  available  for  direct  on-line  measurement. 
Therefore,  the  DAC  is  restricted  to  operate  only  on  information  in  the  available  on-line 
measurements  of  the  system  outputs  and  commands  and  any  disturbance  components  which 
may  actually  be  available  for  direct  measurement.  In  the  case  at  hand,  it  is  assumed  that  none 
of  the  disturbance  components  are  measurable  on-line  and  that  the  information  available 
from  the  plant  consists  of  the  input  command,  command  to  the  actuators  and  the  measured 
pitch  plane  acceleration  and  rate  components  of  the  missile  motion. 


Since  the  idealized  DAC  control  law  is  a  function  of  the  real-time  system  state,  x,  and 
disturbance  state,  z,  the  required  on-line  data  for  practical  DAC  implementation  must  be 
generated  via  use  of  state  reconstructors  (observers)  operating  on  real-time  system  outputs^ 
and  control  inputs  u.  Since  the  external  disturbances  w(t)  are  assumed  to  have  waveform 
structure  and  to  be  modeled  by  known  linear  state  models,  a  state  reconstructor  can  be 
designed  to  generate  estimates  £  of  the  instantaneous  disturbance  state  z.  In  addition,  that 
same  state  reconstructor  can  be  designed  to  produce  estimates^  of  the  instantaneous  system 
state  x. 

Procedures  have  been  developed  to  generate  both  “full-dimensional”  observers  of 
dimension  (n+p),  where  n  is  the  order  of  &  and  p  the  order  of  z,  and  “reduced-dimensional” 
observers  of  dimension  (n+p-m),  where  n,  p  are  as  above  and  m  is  the  rank  of  Q.  The  work 
performed  in  this  study  is  concerned  with  “full-dimensional”  observers. 
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For  the  form  of  the  state  equations  given  by  Equation  ( I ),  the  full-dimensional  observer  is 
expressed  as 


+ 


Jl(t) 


(3) 


where  Kj,  are  gain  matrices  to  be  designed, 

A,  £,  L-  £■  C-  H,  .Q.  M  are  as  previously  described. 


Such  a  composite-type  state  reconstructor  can  be  utilized  to  implement  DAC  control  laws 
in  the  form 
u  =  f(x,z,t). 


Of  course,  for  acceptable  performance  the  real-time  estimation  errors 


4»  =  i4 
tx~lA 


must  settle  to  zero  rapidly  in  comparison  to  system  settling  times  where je,  and  tx  are  given  by 


A  +  FL  +  %(£  +  <&) 

fe  +  %£]  & 

M+KjCC+GL) 

D  +  KjSi 

Slit)/ 


(4) 


3.  PLANT 


The  plant  utilized  for  the  studies  detailed  in  this  report  is  the  pitch  plane  acceleration 
autopilot  channel  shown  in  block  diagram  form  in  Figure  1.  An  ideal  accelerometer  is 
assumed  in  the  acceleration  feedback  loop  and  an  ideal  rate  gyro  is  assumed  in  the  rate 
feedback  loop.  Also,  no  actuator  dynamics  are  considered. 
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Figure  1 .  Pitch  acceleration  autopilot  channel. 


The  transfer  functions  indicated  in  Figure  I  are: 


a.  Pitch  rate  per  fin  deflection  in  pitch: 

0(b)  Kg  (s  +  y) 

—  _ 

6(s)  s  +  2£^  o^s  + 


b.  Lateral  acceleration  per  fin  deflection  in  pitch: 

n(s)  K(s2  +  to  2) 

_  _  n _ n _ 

«s)  s2  +  2ca  +  u^2  * 

with  the  terms  in  the  transfer  functions  being  determined  from  the  aerodynamic 
characteristics  of  the  missile  at  given  points  along  a  trajectory. 


c.  Autopilot  conipensation  terms: 


^(S+a) 
8  +  y 


and 


°n(9**>2 

(8+p)2 


are  compensation  terms 


which  were  designed  into  the  autopilot  to  improve  performance.  Most  of  the  terms  are  varied 
over  a  trajectory  according  to  dynamic  pressure. 
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Since  the  transfer  function  parameters  and  most  of  the  autopilot  compensation  terms  do 
vary  along  the  missile  trajectory,  several  representative  points  along  a  nominal  trajectory 
where  chosen  as  design  points  for  use  in  this  report.  These  points  were  chosen  to  cover  as 
nearly  as  possible  the  entire  range  of  values  of  the  parameters  involved.  Table  /  lists  the  time 
points  and  parameter  values. 

For  the  initial  investigation,  it  was  decided  to  look  at  several  different  configurations 
involving  the  plant,  or  some  part  of  it,  and  a  disturbance  source.  First,  the  entire  loop  was  used 
with  a  disturbance  assumed  to  be  acting  at  the  input.  Next,  the  rate  loop  alone  was  considered 
with  an  assumed  disturbance  summed  in  with  the  b  due  to  fin  deflection.  Third,  the  entire  loop 
was  again  used,  this  time  with  a  disturbance  summed  into  the  output.  As  a  final  case,  the  entire 
loop  was  used  with  disturbances  at  the  input  and  the  output.  Each  of  these  cases  will  be 
detailed  later  in  this  report. 

4.  DISTURBANCE  MODEL 


The  disturbances  modeled  here  in  all  cases  were  taken  to  be  composed  of  constants  plus 
ramps,  i.e.. 


w(t)*  CQ  +  Cxt 


(5) 


where  Co  and  Ci  are,  in  general,  unknown  a  priori  and  can  change  value  in  a  completely 
unknown  random-like  manner.  This  disturbance  model  was  chosen  because  it  is  easy  to  work 
with  but  still  illustrates  the  point. 

To  put  (5)  into  the  form  (2),  proceed  as  follows.  First,  take  the  Laplace  Transform  of  w(t). 


w(s) 


CqS  +  C 


1 


The  characteristic  polynomial  associated  with  this  is 


X 


2 


*  0  . 


(6) 


I! 


O  5  »  8  ftl  ^  " 

©  ai  to  o  cJ  o  s 

i—  i/i  n  w 

CO  CO  o 

■  I  I 


£  8  S 

O  N  f  r  n  8  ' — ’ 

o  ui  8  d  ci  in  s  o 

CM  N  CO  ^  r- 

1  •  •  S 


s  co  n 
8  8  0  8 

cd  d  cd 
T  * 


2  f-  o 

T  ^ 

1  O 


§CO  T 

S  S  N  «  O 

S  i—  Q  W  ft  t- 

d  r-  d  o  in  c\i  t— " 


»  S  c 


d  cm  «p  o  d  ed^  o 

’  •  3 


8  $ 

cm  w  o>  w 

*  8  d  9  8  5 

1  t^3 


5  ° 


<£>  O 

co  ri 

m  oo 

o  P  r>  o 

O  r1  T  ♦  <“ 

"  CO  W  O 

i  i  i  • 


*  *1 


o  N  «■  o 


Therefore,  choose 


w  -  Hz 
S  *  Dz  +  o 


(Note:  no  state  dependence  terms  are  included  in  this  case)  such  that 


JL-  fit 

has  a  characteristic  polynomial  k1  m  0  and  Hi  «as  the  general  form  w  *  Co  +  Cit. 


So,  let 


w-Hz-  (1  0)  (5) 


and 


z  ■  Dz  +  q 


Then, 


dat  Is.  -  Xi | 


"  e2"X  1  I  »  X2  ♦  8  X  «• B,  - 
-e1  -X  1  2  1 


0  .  (7) 


Comparing  (7)  with  (6),  one  must  have  0*  ■  0i  ■  0. 


Thus, 


(i)-[:  :]R) 


♦  a 


w  -  (1 


•ft) 


(8) 
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From  this  one  has,  therefore. 


H  *  (1  0) 


These  two  matrices  are  used  throughout  this  report  in  the  disturbance  state  model. 

In  order  to  give  a  feel  for  just  what  “disturbance”  as  an  entity  is  insofar  as  the  applications 
here  are  concerned,  it  could  be  any  or  all  of:  wind  or  wind  gust,  thrust  misalignment,  tipoff 
rates,  biases,  instrument  drifts,  target  motion  and  more.  An  influencing  agent  which  has 
waveform  structure  and  which  imposes  an  undesirable  effect  on  the  system  may  be  considered 
a  disturbance. 


5.  ACCELERATION  LOOP  WITH  DISTURBANCE  AT  INPUT 
A.  DAC  MODEL  DEVELOPMENT 

A  block  diagram  representation  of  the  autopilot/  disturbance  combination  used  in  the 
development  for  this  section  is  shown  in  Figure  2. 


Figure  2.  Fitch  acceleration  channel  with  disturbance  at  Input. 
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The  closed  loop  transfer  function  for  the  rate  loop  is 


jk> .  *2 

“**>  ?  +  (2Ca^  +  K|C^»  +  to*2  +  Kjc#) 


”0  *  2V*A 

«1  *  2Va  +  K6CR 


*2  *“A  +  kScR 


6(8)  82  +  a08  +  <^2 

-  «  - - : -  . 

Ujfs)  8  +^8  +  82 

With  this,  the  product  of  the  transfer  function  blocks  between  Ui  and  rj  is 
K  Cn{8  +  z) 2 (s2  +w  2)  K  Cn  [s4  +  2*8 3 


(a  +  p)  *  (s*  +  ax8  +  a2)  a’  +  (2p  +  a^s’ 


♦  (z2  +  ta^2)*2  +  2*tan28  +  »n2»2] _ 

+  (p2  +  2pa1  +  a^)*2  +  CajP2  +  2p«2)«  +  a2p' 


b0  -  2z 

2  2 
bl  -  *  ♦  «n 


L 


2  2 
u>  z 
n 


2p  +  a. 


“  P  +  2pax  +  a2 


■  axp  +  2pa2 


The  block  diagram  has  thus  been  reduced  to 


u  +AUl|  K"cnk4  +  V3  +  V2  +  b2s  +  *>3^  * y 


Solving  for  y<s). 


y(s)  -  KnC^  +  i(Kncnbou1(s)  b4y(s)+  J  j Kncnblul(s) 
”b5y(s)  +  5  [KpC^Ujfs)  -  b6y(a) 

+  i  (Kncnb3ul(s)  "  b7y(s*»j)  (12) 


where  j  denotes  an  integration. 

This  can  be  represented  diagrammatically  as 


Figure  3.  Plant  atata  rapraaantatlon. 

and  from  this,  the  equations  for  the  states  can  be  written  directly  as 


*1  •  *2  +  Wo"l  ’  V 

*2  *  *3  *  Kncnblul  ‘  V 
*3  ’  *<  +  Vnb2ul  -  V 


*4  »  KT)^nb3Ul  ”  b7y 


However,  for  purposes  of  DAC  design  these  equations  need  to  be  expressed  as  functions  of  g, 
U.  and  h-  So,  since  Ui  =  u  +  wi. 


U..01M  +  [Kncnl  -  +  tVnl  Si 


Thus,  Equations  (14)  and  (IS)  define  the  remainder  of  the  matrices  needed  for  the  DAC 
design. 

Before  proceeding,  it  is  necessary  to  first  check  for  the  existence  of  a  control,  Uc,  which  can 
totally  counteract  all  the  disturbance  effects. 

This  control  will  exist  if  and  only  if  F  *  B  X  for  someX-  Here, 


K  C 
n  n 


Vb4 

brb5 

b2"b6 

b3"b7 


K  C 
n  n 


b0‘b4 

bl”b5 

Vb6 

b3-b7 


for  r  *  1. 


Such  a  control  does,  therefore,  exist  and  will  be 


H  c  *  “I  Hi  *  ~wi  ”"zi- 


Now,  a  full-dimensional  composite  state  reconstructor  in  the  form  of  Equation  3  must  be 
designed  to  provide  £).  Starting  with  the  error  dynamics  (Equation  4)  in  order  to  obtain  & 
and  fc  we  have  (since  X  *  M  =  0) 


A  +  Kx  C 


(F  +  Kx  G)  H 


D  +  K  2  GH 


To  simplify  notation  in  the  following  development,  let  Equation  (16)  be  represented  as 


e 


A 


ana  let  A  be  represented  as 


eQ  1  0  0 
0  1  0 
e2  0  0  1 

e3  0  0  0 

e4  0  0  0 

e5  0  0  0 


Now,  to  solve  for  the  gain  matrices  ft  and  ft,  one  must  first  find  the  eigenvalues  of  A, 


dot  |A-XI| 

-  C 

ert-X 

1 

0 

0 

e g- 

0 

dot  |A-Xl| 

— 

0 

el 

-X 

1 

0 

0 

®7 

0 

e2 

0 

-X 

1 

e8 

0 

- 

e3 

0 

0 

-X 

e9 

0 

e4 

0 

0 

0 

®10'X 

1 

e5 

0 

0 

0 

®U 

-X 

Expanding  this  determinant  about  the  first  column  results  in  the  expression 


dot  | A- XI |  -  X®- (eQ  +  e1Q)  X5  +  <eoeio"®ll"el“®4®6)  X* 


- 

+  (eOell  +  eIe10-e2-e4e7-«5e6)  X 

+  (ejOj^  +  e2eio“e3-e4e8_e5e7>  X* 

*  (®2®11  +e3®10-e4e9"e5®8)  X 

+  (e3exi~e5®9^  • 

(17) 

21 

9 

■  l 


If  the  desired  roots  of  Equation  (17)  are  Ai.  A2,  Aj,  A*.  A*,  then  the  desired 

characteristic  equation  is 


( X— X x )  (X-X2)  (X-X3)  U-X4)  (A-A5)  (A-Ag)  =  0 


(18) 


Expanding  Equation  (1 8)  and  equating  coefficients  of  like  powers  of  A  between  Equations(l7) 
and  (18)  we  see  that 


(a) 

(b) 

(c) 


e0  +  e10  X1  +  X2  +  x3  +  A4  +  A5  +  Ag 


e0e10“ell“el~e4e6  “  .  E.  .E  .  XiXj  *  A1 

1=1  D=i  +  1 


eOell  +  ele10"e2"e4e7"e5e6  = 


{, 


4  5  6 

*1  I  l  l  XiXiXk 

.i=l  j=i+l  k«j+l  J 


=  -A, 


(d)  eieu  +  e2e10-e3-e4e8-e5e7  » 

3  4  5  6 

£  £  £  Z  A » s 

i=l  j=i+l  k*j+l  l=k+l  1  J  K  1  J 

(e)  e2ell  +  e3e10”e4e9_e5e8  * 

"  I ^1^2^ 3^4 ^  ^  A^) 

^  *3*4*5*$  ^  *2^  ^  *  “A4 

(f)  e3en-e5e9  «  X1X2X3X4X5X6  "  A5  * 
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Substituting  the  relations  for  eo  through  ei  t  from  Equation  (16)  into  (a)  through  (0  and  solving 
for  the  elements  of  Ki  and  K2,  we  obtain 


k 

k 

k 

k 

k 

k 


11 

21 

31 

41 

12 

22 


'KnCnk12  *  b4  *  A0 


"KnCn  <b0k12  +  k22>  +  b5  "  A1 


"KnCn  (b0k22  +  blk12*  *  b«*  A2 


"KriCn  <b2k12  +  blk22’  +  b7  "  A3 


(‘b2KnCnk22  ♦  V  '  Vnb3 


'AS  t  KnCnb3 


09) 


It  is  desirable  that  t(t)  —  0  rapidly,  thus  the  characteristic  roots  Xi,  X2,  Xj,  X4,  Xs,  X«  can  be 
picked  to  best  accomplish  this  depending  on  the  problem  at  hand.  Having  picked  the  X’s,  the 
gain  values  (Equation  19)  can  then  be  calculated.  The  full-dimensional  observer  can  now  be 
implemented,  giving 


*3 

*4 

*1 

*2 


,ku-v 

1 

0 

0 

Wkll+Vb4> 

(k2l“b5) 

0 

1 

0 

Kncn<k2l+brb5> 

(k3rb6> 

0 

0 

1 

knc„‘k3l+b2-b6> 

(k4rb7» 

0 

0 

0 

Knc„<k41+b3-b7» 

k12 

0 

0 

0 

K0Cnk12 

_k22 

0 

0 

0 

KHCnk22 

kl 


kll 

Vn^U+VV 

k21 

KnCn(k21+bl-b5) 

k31 

y  + 

KnCn(k31+b2-b6> 

k41 

KnCn(k41+b3'b7> 

k12 

Wl2 

.  k2i 

*ncnk22 

(20) 
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Figure  4  is  a  diagram  of  the  composite  plant-DAC  system.  On  the  diagram,  hi  •  hj  are  the 
components  of  the  last  matrix  on  the  right-hand  side  of  (Equation  20).  The  remainder  of  the 
symbols  have  been  previously  defined. 

So  now  we  have  a  disturbance  cancelling  control  term,  Uc  =  -zi,  and  we  have  a  composite 
state  reconstructor  which  gives  zi.  The  questions  now  are: 

•  Can  the  X’s  be  picked  so  that  t,  and  t,  settle  to  zero  “rapidly”? 

•  If  so,  does  Uc  really  cancel  out  the  effects  of  wi?  If  both  of  these  can  be  answered  in  the 
affirmative,  then 

— How  well  does  the  DAC  work  if  the  plant  parameters  are  varied  from  the  design 
point? 

— How  do  the  DAC  characteristics  vary  over  the  trajectory? 

Simulation  results  should  provide  answers  to  these  questions. 

B.  SIMULATION  AND  RESULTS 

The  composite  system  shown  in  Figure  4  was  simulated  on  a  digital  computer.  The 
simulation  was  written  so  that  the  plant  parameters  could  be  arbitrarily  varied  around  the 
point  for  which  the  DAC  was  designed.  A  listing  of  this  simulation  is  given  in  Appendix  A. 

As  a  first  cut  at  seeing  how  effective  the  DAC  would  be,  several  runs  were  made  for  the  t  = 
9.85  sec  and  t  =  18  sec  points.  Figures  5  and  6  give  the  results.  As  can  be  seen,  the  DAC 
effectively  cancels  out  an  input  disturbance  (wi  =  1.0)  equal  to  the  input  command. 

To  check  the  sensitivity  of  the  DAC  to  plant  parameter  variations,  a  series  of  runs  were 
made  with  parameters  varied  around  the  t  =  1 8  sec  values.  Table  2  is  a  summary  of  the  results 
obtained  and  Figures  7  through  40  give  the  system  output,  y,  and  reconstructor  state,  zi 
(disturbance  estimate),  for  each  case.  The  table  shows  how  the  peak  value  of  y  varied  and  how 
the  peak  value  and  settling  time  of  zi  varied  due  to  both  individual  and  collective  parameter 
changes.  In  all  cases,  the  input  command  is  I  and  the  settling  time  is  defined  to  be  the  time  at 
which  the  response  stays  within  ±5%  of  steady-state. 
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Figure  9.  System  output  response  (y),  W-|  =  1,  +10S  vi 


Figure  12.  OAC  disturbance  estimate^),  W-j  =  1,  -10%  variation  on  K£ 
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=  *-M  ‘(A)  MuodMJ  indjno  UMftAs  *61-  AinBu 


Figure  20.  DAC  disturbance  estimate  (z-|),  W-*  =  1,  -10%  variation  on 


Figure  21.  System  output  response  (y),  Wi  =  1,  +10%  variation  on 


Figure  23.  System  output  response  (y),  W-|  =  1,  -10S  variation  on  »*. 
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Figure  29.  System  output  response  (y),  W-j  =  1,  +10%  variation  on  Cr. 


System  output  response  (y),  Wj  =  1,  -10%  variation  on  Cr. 


Figure  96.  DAC  disturbance  estimate  (z-f),  Wf  = 
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System  output  response  (y),  Wf  =  1,  +20%  variation  on  aN 


Figure  38.  DAC  disturbance  estimate^),  Wi  =  1,  +20%  variation  on  all 
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Figure  39.  System  output  response  (y),  Wf  =  1,  -20%  variation  on  aH 
parameters. 


Figure  40.  DAC  disturbance  Mtfmata(f i),  Wj  = 
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From  these  results,  the  first  two  questions  previously  posed  can  indeed  be  answered  in  the 
affirmative.  The  answer  to  the  third  question  would  seem  to  be  that  the  DAC  works  well  with 
at  least  up  to  20%  variation  of  plant  parameters  and  possibly  for  larger  variations.  For  a  given 
system,  though,  this  should  be  thoroughly  verified  by  checking  at  all  critical  times  along  a 
trajectory,  i.e.,  burnout,  apogee,  etc. 


In  order  to  answer  the  fourth  question,  three  of  the  time  points  shown  in  Table  /  were  used 
in  the  simulation.  The  roots  of  Equation  (18),  which  were  used  to  settle  out  the  state 


TABLE  3.  ROOTS  FOR  DETERMINING  DAC  GAIN  MATRICES 


reconstructor  in  each  case  and  to  calculate  the  components  of  the  DAC  gain  matrices,  are 
shown  in  Table  3.  The  components  of  &  and  &  are  shown  in  Table  4. 


TABLE  4.  DAC  GAIN  MATRIX  COMPONENTS 


Three  simulation  runs  were  made  at  each  time  point: 

•  with  nominal  airframe  parameters,  no  disturbance, 

•  with  nominal  airframe  parameters  and  a  disturbance  and 

•  with  a  20%  variation  on  airframe  parameters  in  the  direction  of  increasing  flight  time, 
with  a  disturbance.  The  results  are  presented  in  Figures  41  through  58. 

From  these  results  and  the  DAC  parameters  shown  in  Tables  3  and  4,  it  is  evident  that  a 
DAC  designed  at  one  point  of  a  trajectory  will  not  perform  as  well  as  needed  over  large 
portions  of  the  trajectory.  Gain  switching,  similar  to  an  autopilot  gain  switch  program,  wiV.  be 
required  for  DAC  implementation. 

C.  CONCLUSIONS 

For  this  case,  with  the  disturbance  at  the  input,  it  was  possible  to  find  a  control  Uj  which 
could  be  implemented  and  which,  theoretically,  would  totally  cancel  the  disturbance.  In  a 
practical  application,  it  was  found  that  the  control  did  cancel  the  disturbance  very  well,  that 
the  DAC  would  continue  to  function  well  within  a  band  about  the  design  point  and  that,  with 
gain  switching,  the  DAC  should  perform  its  function  as  the  plant  parameters  vary  over  an 
entire  trajectory.  As  can  be  seen  from  Table  4  the  DAC  gains  do  have  a  wide  range. 

For  the  apogee  case,  since  the  system  is  so  sluggish,  the  DAC  did  not  offer  much  in  the  way 
of  disturbance  cancellation,  i.e.,  the  estimation  errors  did  not  settle  out  rapidly  enough.  One 
reason  is  due  to  the  nearness  of  the  eigenvalues  of  the  A,  matrix  to  zero.  This  allowed  the  large 
overshoot.  However,  these  eigenvalues  had  to  be  maintained  in  this  region  because  moving 
them  to  more  negative  positions  caused  an  instability  to  develop.  Overall,  it  might  be 
advantageous  to  zero  out  the  DAC  gains  near  apogee. 

6.  RATE  LOOP  WITH  DISTURBANCE  ON  OUTPUT 
A.  DAC  MODEL  DEVELOPMENT 

The  missile  from  which  this  autopilot  channel  was  taken  uses  an  attitude  control  during 
boost,  so  it  was  of  interest  to  consider  the  rate  loop  alone,  with  a  disturbance  included,  to  see  if 
a  DAC  would  be  useful  in  taking  out  effects  due  to  externa)  rate  perturbations. 
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Figure  43.  System  output  reeponee,  tf  =  9.85  sec,  POO  =  1,  W-j  =  1.0. 


Figure  44.  OAC  disturbance  estimation  error,  tf  =  9.85  sec,  PQO  =  1,  W1  =  1.0. 


Figure  45.  System  output  response,  tf  =  9.85  eec,  W-|  =  1.0, 
airframe  parameters. 


Figure  46.  DAC  disturbance  estimation  error,  tf  =  9.85  sec,  Wj  -  1.0, 
-20%  variation  on  airframe  parameters. 


response,  tf  =  66.7  tec,  PGO  =  0.5,  Wf 


Figure  49.  System  output  response,  tf  =  68.7  sec,  PGO  =  0.5,  W-j  =  0.5. 


Figure  50.  DAC  disturbance  estimation  error,  tf  =  66.7  sec, 
Wi  =  0.5. 


i  output  response,  tf  =  66.7  esc,  FQO  =  0.5,  Wi  =  0.5, 
variation  on  airframe  parameters. 


Figure  52.  DAC  disturbance  estimation  error,  tf  =  65.7  sec,  PGO  =  0.5, 
Wi  =  0.5,  +.20%  variation  on  airframe  parameters. 
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1.0,  +20%  variation  on  airframe  parameters. 


The  rate  loop  shown  in  Figure  2  was  rearranged  and  simplified  as  shown  in  Figure 59.  Thi 
compensation  term  was  reduced  to  Cr  since  no  actuator  dynamics  are  considered.  Th 
disturbance  is  shown  as  a  rate  imposed  on  the  airframe  in  addition  to  that  due  to  a  given  fli 
deflection.  Therefore,  the  total  body  rate. 


9m  “  ®A/F  +  w2  . 

From  Figure  59  one  has 

§A/F(S)  _  K0(8+1r) 

ss  -  :  5 

6(s)  s  +2c^u^S+o)^  • 

Therefore, 

s2eA/F<s)  +  2tA»A.9Vp(»)  +  .j|  5VP(S)  - 

K»s6 (b)  +  Kgy6 (s) 
and 


0A/F(s) 


?  (Kii(,)-2!j“A/P,s,t  S  (kjt« (s)-ba6a/f(,)1 
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This  appears,  along  with  the  rest  of  the  loop,  as 


Writing  the  state  space  equations  directly  from  this  and  substituting  for  £  gives: 

X2  =  K0y6“  “A  X1  =  -(K6YCR  "  “A)  X1  +  K6yUR  “  KeYCRW2 

X1  =  x2  +  Kg6-  “  (ReCR  +  25A“A)  X1 

+  x2  +  KgUR  -  KgCRw2 

Y  ■  (*1  +  w2)  CR 

x  (2 

6  =  uR  -  y  . 

Expressing  these  in  the  form  (1), 


Z"  [CR  0]  +  [0]UR+  [CRJ  w2 


Now;  in  this  case,  £  =  £  £  for  £  =  [-Cr  ).  Therefore,  the  theoretical  total  absorption  controller 
would  be  Uc  =  |Cn]  ^2.  If  this  is  implemented,  however,  it  does  not  remove  w2  from  the  output, 
y[.  But,  the  desired  action  in  this  case  is  to  remove  any  disturbance  rates  imposed  on  the  missile 
so  that  it  will  maintain  a  given  attitude  during  boost.  On  examining  the  plant,  it  can  be  seen 
that  the  disturbance  rate  %2  imposed  on  the  body  can  be  related  to  an  “equivalent**  fin 
deflection  6d,  i.e„  it  can  be  considered  as  an  additional  body  rate  which  would  have  resulted 
from  an  additional  fin  deflection  command  fiD.  In  other  words, 

=  eA/F  +  w2  =  (6  +  V  I  * 

So,  if  the  proper  gain  can  be  found,  and  if  a  state  observer  can  be  designed  which  will 
reconstruct  zi,  then  Uc  =  -CRGzi  can  perhaps  be  used  as  a  partial  absorption  (minimization) 
control  on  the  effects  of  »2,  where  Crg  is  the  sought  gain. 


Again,  Equation  (4)  is  used  to  obtain  the  gain  matrices  Ki,  In  this  case, 

“  + 2  c^oi.)  1  k,  ,  1  I 

6  R  A  A  +  11  [CR  O]  11  0  |  [CR  OJ  | 

;(Vcr+“a)  °J  iNJ  k2i-V 


kij 

ICI 

k  1 

_22_ 


o  i]  h12 

+  11  [c. 

_°  °J  Lk22j 


e+/2 


;llCR_2cA“A"KeCR  1 

(kn-1 

21CR'"A'KeYCR  0 

(k2l"' 

12CR  0 

k12CR 

22CR  0 

k22CR 
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Let  Equation  (23)  be  written 


Be  +  *  o' 


and  g  be  written  as 


B  * 


B, 


B, 


B- 


B , 


1 

0 

0 

0 


B5  0 

B6  0 
B7  1 
B8  0 


Solve  for  the  eigenvalues  of  |E 
det  |B  -  Al|  -  0 


det  |B  -  Al|  - 


Bj-A 

Bo 


B- 


B> 


1 

-A 

0 

0 


Br 


B< 


B?-A 

B0 


0 

0 

1 

-A 


*  0 


This  can  be  expanded  to  give 


a4  -  (b1+b3) a3+ (b1b3-b3b4-b3b5-b2)  A2 


+  (B^-BjBg+BjBj-B^g)  A 


+  (b4b2-b6b4)  ■  0  . 


(24) 


If  the  desired  roots  of  Equation  (24)  are  Ai,  X2,  X3,  X4,  then  the  desired  characteristic  equation  is 


<A-A1)  (A-A2)  (A-A3)  (A-A4)  =  0 


(25) 
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Expanding  Equation  (25).  equating  coefficients  of  like  powers  of  A  between  Equations  (24) 
and  (25).  substituting  back  in  for  Bi  through  B»  and  solving  for  ku.  k2),  ki2,  k22  gives 


kll  *  (25AwA  “  k12CR  +  Xl+X2+X3+X4  +  KeCR)/,<:R 


k21  "  I2?A“ACR(k12“k22)  +  kllk12CR  ”*A  "K6k22CR 


+  X1X2  +  X1X3  +  X1X4  +  X2X3  +  X2X4  +  X3X4^  ^  (“C^) 


k12  "  (“2tA“ACRk22  +  X1X3X4  +  X2X3X4  +  X1X2X3 


+  x1x2x4)  /  (C8«“) 


22 


X1X2X3X4  ^  CRWA 


Again,  the  X’s  are  picked  such  that  e(t)  —  0  rapidly.  With  these,  the  gains  can  be  determined. 
The  full-dimensional  observer,  in  the  form  of  Equation  (3),  for  this  case  is 


*1 

kllCR  “ 

2ca“a 

1 

kllCR 

o" 

*1 

kn 

V 

*2 

k21CR  “ 

2 

“A 

0 

k21CR 

0 

^2 

k21 

y  + 

V 

*1 

k12CR 

0 

k12CR 

1 

A 

*1 

k12 

o 

k 

k22CR 

0 

k22CR 

0 

L2*J 

_k2i 

0 

(26) 


So,  the  question  posed  here  is:  Can  the  proposed  gain  Cro  be  found  so  that,  in  conjunction 
with  the  state  reconstructor,  the  effects  of  Ha  can  be  minimized? 

For  answers  to  these  questions,  a  simulation  is  again  used. 
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B.  SIMULATION  AND  RESULTS 


Figure  60  is  a  diagram  of  the  composite  plant-D  AC  system  which  was  simulated  on  a  digital 
computer.  A  listing  of  this  simulation  is  shown  in  Appendix  B. 

The  controller  in  this  loop,  with  disturbance  as  shown,  would  probably  be  used  only  during 
the  attitude  controlled  boost  phase  of  this  missile.  However,  for  illustrative  purposes,  two  of 
the  time  points  shown  in  Table  I  (9.85  sec,  135.8  sec)  will  be  used  for  this  investigation. 


Figures  61  through  72  present  the  results  obtained  for  the  time  points  above.  For  each  point, 
three  runs  are  presented:  (1)  nominal  run  with  no  disturbance;  (2)  run  with  W2  equal  to  the 
plant  steady  state  response  (x()  due  to  the  input  command,  from  (1),  but  with  CRQ  =  0.;  (3) 
same  as  (2)  except  Crg  is  given  an  appropriate  value.  This  value  is  determined  from  the  ratio 
(5/  xi)»  from  the  undisturbed  case,  as  would  be  expected. 

For  the  9.85  sec  point,  Crg  -  -3.54,  and  it  can  be  seen  from  comparison  of  Figures  63  and  65 
that  the  effects  of  w2  are  largely  removed.  For  the  t  *  135.8  sec  point,  Crg  =  -1.45,  and  similar 
conclusions  are  reached  (compare  Figures  69  and  71). 

The  gain  matrix  components  for  the  DAC  in  the  two  cases  above  are  given  in  Table5vthen 
the  eigenvalues  of  matrix  B  were  taken  to  be  Xi,  =  -3.,  X2  =  -5.,  Xj  =  -4.  +  jl.,  X4  =  -4.- jl. 

TABLE  5.  DAC  GAIN  COMPONENTS 
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AD-A083  765 

UNCLASSIFIED 


ARMY  MISSILE  COMMAND  REDSTONE  ARSENAL  AL  TECHNOLOGY  LAB  F/6  17/7 
INVESTIGATION  OF  DISTURBANCE  ACCOMMODATING  CONTROLLER  APPLICATI— ETC f U) 
MAY  79  N  L  HCCOKAN 

DRSM I -T-79-63  NL 


MICROCOPY  RESOLUTION  T t ST  CHART 

NAIIONAl  HURL  All  01  S1AN0AK0S  l«»M  A 


TIME  (sec) 

Figure  66.  OAC  disturbance  estimation  error,  tf  -  9.65  sec,  Ur  =  1,  Wg  = 
•0.26,  Crq  =  -3.54. 


•ix 


•  «  TIME  (mc)  •  ■ 
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TIME  (tec) 

Figure  69.  Rate  loop  reeponse,  tf  =  1 35.8  sec,  Ur  =  1 ,  W2  =  -0.5588,  Crq  = 


TIME  (mc) 

OAC  dteturbanc*  — MwHon  wror,  If  =  13SJ  — c,  Ur 
-0.5588,  Crq  =  0. 


Figure  71.  Bate  loop  response,  tf  =  135.8  sec,  Ur  =  1,  Wg  =  -0.5588,  Crqc-1.< 


Figure  72.  DAC  disturbance  estimation  error,  tf  =  195.8  sec,  Ur  =  1,  Wj  s 
-0.5588,  Con  =  -1.45. 


C.  CONCLUSIONS 


From  the  results  obtained  here,  even  though  there  was  no  total  absorption  control  u*  which 
would  exactly  cancel  »2,  by  using  a  state  reconstructor  to  estimate  the  value  of  the 
disturbance  an  implementation  was  possible  whereby  the  effects  of  w2  could  be  minimized. 

7.  ACCELERATION  LOOP  WITH  DISTURBANCE  ON  OUTPUT 

A.  DAC  MODEL  DEVELOPMENT 

As  with  the  rate  loop  development  in  Section  6,  this  section  will  consider  again  a 
disturbance  summed  with  the  output,  this  time  for  the  acceleration  loop.  This  case  is  of  much 
interest  since  the  acceleration  autopilot  is  used  in  the  control  loop  from  burnout  onwards.  A 
block  diagram  of  this  loop  is  shown  in  Figure  73. 


Ur  ±r\  u 

KnCn [s4+b0s3+b1s2+b2s+b3 

"2+ftY 

=  nm 

- 

j 

- 2 - g - 5 - 

s  +b4s  +b^s  +bgS+b7 

Figure  73.  Acceleration  loop  with  disturbance  on  output. 


The  transfer  function  from  u  to  u2  can  be  represented  identically  as  shown  in  Figure  3  with 
the  same  states  and  parameters.  However,  in  this  case,  the  matrix  representation  in  the  form  of 
Equation  (1)  will  be  written  as 


-b4  1 


-b5  0 


"b6  0 


-b?  0 


0 

1 

0 

0 


X1 

'J* 

XX 

i 

o 

X2 

+  K  C 

bl-b5 

n  n 

X3 

b2"b6 

LX4. 

b3-b7. 

u  +  [0]  w3 

(27) 
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l  i  o  o  o  ]  ♦  lKncni  u  +  m  w3 


Using  the  same  approach  in  this  case  as  was  used  with  the  rate  loop  in  the  previous  section,  and 
for  the  same  reason,  a  gain  (Cral)  is  sought  which  can  be  used  in  conjunction  with  a 
disturbance  state  reconstructor  output  to  accomplish  a  minimization  of  the  effects  of  on  y. 


Proceeding  as  before  with  the  £  dynamics  one  has 


-b4  l  o  o 
-b5  0  1  0 

-b6  0  0  1 

-b7  0  0  0 


[100  0] 


[1]  [1  0] 


o 

£  +  O 


k12 

U'1 


0  0  0] 


I.  -I 

0  1  It- 2 

+  12  [1  0] 

1°  M 


which  reduces  to 


(kll“b4) 

1 

0 

0 

kll 

(Jt21-b5) 

0 

1 

0 

k21 

<k3l"b6) 

0 

0 

1 

k31 

(k41~b7) 

0 

0 

0 

k41 

k12 

0 

0 

0 

k12 

k22 

0 

0 

0 

k22 

-  C  e  +  ^ 
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Solving  for  the  eigenvalues  of  ?  gives 


X6  +  *  +  ^b5~b4k12~k22~k21)  ^ 

+  (b6-b4k22-b5k12_k31)  X  +  (b7"b5k22~b6k12~k41)  ^ 

+  ^-b6k22-b7k12^  *  ~  b7k22  =  0  (29) 

In  the  same  manner  as  for  the  first  two  cases,  one  can  solve  for  the  components  of  the  gain 


matrices,  Ki 

and  K2.  Doing  so  gives  the  following, 

(a) 

kll  = 

b4-k12  +  A0 

(b) 

k21  = 

b5-b4k12"k22-Al 

(c) 

k31  = 

b6-b4k22~b5k12  +  A2 

(d) 

ii 

rH 

b7'b5k22‘b6k12"A3 

(e) 

k12  = 

<b6k22-V  /b7 

(f) 

k22  = 

-A5/b7 

where  At  through  Aj  are  as  defined  for  use  in  Equation  ( 19).  The  full  dimensional  observer  can 
now  be  expressed  as 
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The  question  for  this  case  is  of  the  same  type  as  for  the  rate  loop,  i.e.,  can  a  gain,  Crai.,  be 
1  found  and  used  inconjunction  with  the  slate  reconstructor  output  zi,  to  minimize  the  effects  of 
the  disturbance? 


B.  SIMULATION  AND  RESULTS 

The  diagram  for  this  composite  system,  with  proposed  control,  is  shown  in  Figure  74.  The 
e’s  and  h’s  are  as  defined  for  Figure  3  with  the  following  exceptions: 


e6  " 

kll 

A 

■~4 

II 

k21 

e8  ~ 

k31 

ll 

<T\ 

0 

k41 

A 

I-1 

O 

II 

'  k12 

G11  " 

:  k22 

A  listing  of  the  digital  simulation  is  given  in  Appendix  C. 

For  this  loop,  the  gain,  Cral,  is  determined  initially  from  the  ratio  y  (see  Figure  74)  from 
the  undisturbed  case  and  is  then  iterated,  if  necessary,  to  obtain  a  final  value. 


Several  of  the  time  points  from  Table  I  were  used  to  analyze  this  case.  Figures  75  through  84 
give  results  for  the  9.85  sec  airframe  parameters.  By  comparing  Figures  75,  77  and  79,  it  can  be 
seen  that  the  effects  of  the  disturbance  are  cancelled  for  a  disturbance  magnitude  equal  to  the 
input  command.  Figures  81  and  82  show  results  for  wj  equal  to  twice  the  input  command 
magnitude  and  Figures  83  and  84  are  results  for  wj  equal  to  five  times  the  input  command. 

For  tr  =  66.7  sec  (apogee),  an  input  command  of  0.5  was  used  with  W3  -  0.5.  Figures  85 
through  89  give  the  results  for  this  time  point.  As  can  be  seen  from  Figures  85, 87  and  89,  the 
disturbance  effects  were  again  cancelled  out  although,  since  the  system  is  sluggish,  it  takes 
longer  to  settle  out  than  the  previous  case. 
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Acceleration  loop  reeponee,  tf  =  9.85  sec,  PGO  =  1,  W3  = 


■  *  TIME  (mc) 

Fig ura  77.  Acceleration  loop  raaponM,  tf  =  9.85  mc,  PGO  =  1,  W3 
CRAL  =  0. 


Figure  78.  DAC  disturbance  estimation  error,  tf  =  9.85  sec,  PGO  =  1,  W3  = 
1.  CraL  = 


Figure  79.  Acceleration  loop  reeponee,  tf  =  9.85  tec,  POO  =  1,  W3  = 
Oral  =  -0.12. 


Figure  P5.  Acceleration  loop  reeponse,  tf  =  66.7  eec,  POO  =  0.5,  W3  = 


W3  =  0.5,  CraL  = 


Figure  89.  Acceleration  loop  reeponae,  ty  =  66.7  sec,  PGO  =  0.5,  W3 = 0.5, 
CRAL  =  -1-50. 


m 


For  ti  =  1 1 1.4  see.  an  input  command  and  disturbance  magnitude  of  1.0  were  again  used. 
Figures  90  through  94  show  the  results  obtained.  Comparing  Figures  90.  92  and  94,  it  can  be 
seen  that  the  distrubance  cllecis  are  removed.  The  output  on  Figure  94  could  perhaps  be 
settled  out  better  by  more  iterations  with  the  state  reconstructor  roots. 

Table  6  gives  the  components  of  Kiand  K2and  the  roots  of  Equation  (29)  used  in  the  above 
three  cases  along  with  the  value  for  Crai.  in  each  case. 


TABLE  6.  STATE  RECONSTRUCTOR  DATA  AND  CRAL  FOR  ACCELERATION 
LOOP  WITH  DISTURBANCE  ON  OUTPUT 


\  TIME 

POINT 
\^SEC) 

PARAMETER'S^ 

9.85 

66.7 

111.4 

k11 

3.653 

-0.129 

-14.54 

k21 

195.44 

-2.033 

-184.34 

k31 

2159.92 

-3.68 

-605.05 

k41 

2668.13 

-1.776 

-483.25 

k12 

-12.26 

-4.584 

7.80 

k22 

0.0 

-2.793 

-29.47 

X1 

-3. 

-1. 

-2. 

x2 

-3. 

-1. 

-2. 

X3 

-4  +  j4 

-1.5  +  j0.25 

-3  +  jl 

x4 

-4  -  j4 

-1.5-j0.25 

-3  -  jl 

X5 

-6. 

-2. 

-5. 

X6 

-6. 

-2. 

-5. 

Oral 

-0.12 

-1.50 

-0.40 
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Figure  90.  Acceleration  loop  reeponee,  tf  =  111.4  see,  PGO  =  1,  Wg 
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Figure  92.  Acceleration  loop  response,  tf  =  111.4  sec,  PGO  =  1,  W3  = 
CRAL  = 


Figure  94.  Acceleration  loop  response,  tf  =  111.4  sec,  PGO  =  1,  W3  = 

CraL  =  -°-4- 


C.  CONCLUSIONS 


From  the  results,  in  similar  fashion  to  the  rate  loop,  even  though  no  u*  existed  which  would 
exactly  cancel  wj,  by  using  a  state  reconstructor  to  estimate  the  value  of  the  disturbance  an  im¬ 
plementation  was  possible  whereby  the  effects  of  Wj  could  be  minimized.  Again,  gain  switch¬ 
ing  would  be  required  to  implement  this  in  a  system. 

8.  ACCELERATION  LOOP  WITH  INPUT  AND  OUTPUT 
DISTURBANCES 

A.  MODEL 

As  a  last  case  for  this  report,  the  acceleration  loop  with  m  and  1*3  both  included  is 
considered.  If  the  procedures  given  in  previous  sections  are  followed  in  attempting  to  derive  a 
D  AC  for  this  case,  it  becomes  necessary  to  evaluate  the  determinant  of  an  8  X  8  matrix  to  solve 
for  the  components  of  the  gain  matrices  Ki  and  K:.  This  evaluation  is  tedious  at  best  with 
many  opportunities  for  mistake. 


Therefore,  it  is  of  interest  to  see  if  the  DACs  developed  in  Sections  5  and  7  can  be  combined 
in  such  a  fashion  as  to  continue  to  function  as  desired  in  cancelling  the  effects  of  m  and  A 
block  diagram  of  the  proposed  combination  is  shown  in  Figure  95. 


Figure  95.  Qeneral  block  diagram  of  plant/DAC  with  both  acceleration 
dleturbance  Inputs. 
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In  this  development,  the  same  plant  state  equations  as  before  are  used  with  changes  only  in 
nomenclature  where  necessary.  The  basic  state  reconstructor  models  are  as  previously 
developed.  In  this  case,  however,  the  rearrangement  of  the  plant  output  portion  of  the  data 
input  to  the  wi  state  reconstructor  should  be  noted.  Since  this  DAC  was  developed  with  no 
disturbance  on  the  plant  output,  in  order  for  it  to  function  properly  it  is  necessary  to  use  a 
plant  output  with  Wj  removed.  This  is  possible  since  the  wj  state  reconstructor  is  also 
reconstructing  the  plant  states.  So,  where  in  Section  5  the  plant  output  is, 

y  "  Kncnul  +  X1  • 
here  it  is  formulated  as 

*PR  *  Kncnul  +  *13  • 

Thus,  it  is  important  in  this  application  for  the  Wj  reconstructor  to  settle  out  as  rapidly  as 
possible. 

B.  SIMULATION  AND  RESULTS 

A  listing  of  the  simulation  is  given  in  Appendix  D.  In  this  particular  case,  only  one  time 
point  from  Table  I  (9.8S  sec)  was  used  since  the  purpose  here  was  to  see  if  two  DAC’s  could  be 
operated  successfully  in  a  serially  connected  mode. 

Figures  96  through  98  show  the  loop  output,  y,  and  the  disturbance  estimation  errors, 
and  e<^,,  for  a  nominal  run,  i.e.,  input  of  1.0,  no  disturbances.  This  output  compares  with 
similar  case  results  from  Section  5  as  would  be  expected.  In  order  to  check  each  reconstructor 
and  see  if  any  undesirable  interactions  were  taking  place,  two  runs  were  made,  one  with  wi  = 
1  „  W3  =  0.  and  one  with  Wi  -  0.,  W]  =  1.  The  results  are  shown  in  Figures  99  through  104.  In 
the  first  case,  everything  looks  okay.  In  the  second  case,  since  the  W3  reconstructor  is  feeding 
input  to  the  wi  reconstructor  and  has  a  settling  time  of  several  seconds,  there  are  some 
dynamics  induced  in  the  wi  reconstructor.  This  in  turn  causes  some  dynamics  to  appear  in  the 
output.  However,  if  Figure  102  is  compared  to  Figures  96  and  99  on  a  similar  scale,  the  results 
do  not  have  such  an  undesirable  appearance.  From  this  it  can  be  seen  that  some  interaction  is 
taking  place  but  not  to  such  an  extent  that  the  DAC  performance  in  either  case  is  impaired. 

The  remainder  of  the  runs  were  made  with  disturbence  inputs  on  w^  and  wj 
simultaneously.  Figures  103  through  107  give  results  for  wi  *  1.,  w»  *  2.;  Figures  108  through 
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HO  for  wi  =  I.  +  0.2t,  wi  =  0.5  +  0. It  and  Figures  III  through  113  for  wt  =  l.-0.2t,  w3  =  0.5 
+  0.5t.  In  all  these  cases,  even  though  the  induced  dynamics  are  noted  in  «**,,  the  DAC’s  per¬ 
formed  their  function  of  cancelling  the  disturbance  effects. 

C.  CONCLUSIONS 

From  the  results  in  this  section  it  would  appear  that  it  is  possible  to  design  DAC’s  for 
separate  disturbances  in  different  parts  of  a  loop,  thereby  simplifying  the  size  of  the  matrices 
involved  in  the  calculations,  and  combine  them  in  a  simple  manner  to  achieve  the  desired 
results. 

9.  CONCLUSIONS 

Conclusions  have  been  presented  in  Sections  5  through  8  regarding  the  results  obtained  with 
the  design  in  each  section.  Overall,  it  has  been  shown  that  it  was  possible  to  cancel  out  or 
minimize  the  effects  of  the  disturbances  modeled  herein  by  use  of  DAC  techniques.  It  was  also 
shown  that  it  was  possible  to  combine  two  separate  DAC's,  designed  for  disturbances  at 
different  places  in  the  plant,  into  a  functioning  unit  which  would  still  perform  its  overall 
purpose.  This  is  especially  important  since  the  size  of  the  matrices  involved  in  designing  a  “full- 
dimensional”  observer  is  directly  related  to  the  dimensions  of  the  plant  plus  disturbance 
models.  Thus,  any  procedure  which  can  reduce  the  dimensionality  involved  is  important.  In 
this  regard,  several  of  the  designs  here  might  be  redone  utilizing  a  “reduced-order"  observer  to 
see  how  well  such  a  DAC  would  perform. 

Although  it  would  appear  from  the  results  obtained  here  that  a  DAC  might  be  very  useful  in 
cancelling  out  unwanted  disturbances,  the  only  way  to  really  be  sure  how  one  would  function 
in  a  system  application  would  be  to  implement  one  in  a  6-DOF  simulation  and  fly  it  with  a 
severe  program  of  varying  disturbance  vectors. 
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Figure  96.  Acceleration  loop  response.  If  =  9.85  eec,  PGO  =  1,W-j  =  Wg  = 


XI* 


•  *  4  TIME  (sec)  • 

Figure  98.  OAC  disturbance  estimation  error  for  Wj,  tf  =  9.85  sec,  PGO  = 


Figure  96.  Acceleration  loop  response  tf  =  9.85  sec,  PGO  =  1,  W-j  =  1,  W3  = 


Figure  100.  DAC  disturbance  estimation  error  for  W-j,  If  =  9.85  NO,  POO  =  1,  W-|  *  1,  W3 


1.0 ,  W3  = 


Figure  102.  Acceleration  loop  response,  tf  =  9.85  sec,  POO  =  1,  W-j  =  0,  W3 


Figure  106.  DAC  disturbance  estimation  error  for  Wi,  tf  =  9.85  sec,  PQO  =  1,  W-j  * 
W3  =  2. 


W3  = 


TIME  (mc) 

Figure  108.  Acceleration  loop  reeponoe,  tf = 9.85  eec,  PGO  =  1,  Wj  =  1  +  0.2t,  W3  =  0.5 

+  0.1t. 


Figure  110.  DAC  disturbance  estimation  error  In  W3,  tf  =  9.85  sec,  PGO=  1,  W-j  =  1  + 
0.21,  W3  =  0.5  +  0.11. 


TIME  (sec) 

Figure  112.  DAC  disturbance  estimation  error  for  W-| ,  tf  =  9.85  sec,  PGO-1,W<|  =  1.0 
-  0.2t,  W3  =  0.5  +  0.5t. 


Figure  113.  DAC  disturbance  estimation  error  for  W3,  tf = 
-  0.2t,  W3  =  0.5  +  0.5L 
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APPENDIX  A 

DIGITAL  SIMULATION  OF  ACCELERATION 
LOOP  WITH  DISTURBANCE  ON  INPUT 
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APPENDIX  B 

DIGITAL  SIMULATION  OF  RATE  LOOP 
WITH  DISTURBANCE  ON  OUTPUT 
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APPENDIX  C 

DIGITAL  SIMULATION  OF  ACCELERATION 
LOOP  WITH  DISTURBANCE  ON  OUTPUT 
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APPENDIX  D 

DIGITAL  SIMULATION  OF  COMPOSITE  ACCELERATION 
LOOP  WITH  DISTURBANCE  ON  INPUT  AND  OUTPUT 
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